Cosmological Parameters from the SDSS DR5 Velocity Dispersion 
Function of Early-type Galaxies through Radio- Selected Lens 
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ABSTRACT 

We improve strong lensing constraints on cosmological parameters in light of 
the new measurement of the velocity dispersion function of early-type galaxies 
based on the SDSS DR5 data and recent semi-analytical modeling of galaxy for- 
mation. Using both the number statistics of the CLASS statistical sample and 
the image separation distribution of the CLASS and the PANELS radio-selected 
lenses, we find the cosmological matter density f2 m>0 = 0. 25^0^8 (68% CL) as- 
suming evolutions of galaxies predicted by a semi-analytical model of galaxy 
formation and fi m ,o = 0.26jZo!o8 assuming no evolution of galaxies for a flat cos- 
mology with an Einstein cosmological constant. For a flat cosmology with a 
generalized dark energy, we find the non-evolving dark energy equation of state 
w x < -1.2 (w x < -0.5) at the 68% CL (95% CL). 



Subject headings: cosmological parameters 
statistics — gravitational lensing 



galaxies: evolution — galaxies: 



1. Introduction 

S trong lensing has been an important astrophysical too l for probing bot h cosmology 
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is also potentially a useful tool to test theories of gravity (e.g., 
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2005, 


2006) 
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At the time of this writing there are ~ 90 galactic-scale strong lensesQ Parts of them 
form well-defined samples that are useful for stati stical analyses. For example, 26 lenses 



from the Cosmic Lens ALL-Sky Survey (CLASS; iMvers et al. 



2003 



. 3rowne et al.l 120031 ) 



Winn et al 



200l|) form a well- 



and the PMN-NVSS Extragalactic Lens Survey (PANELS; 
defin ed radio-selected lens sample, and the Sloan Digital Sky Survey (SDSS; IQguri et al. 



20061 ) has accumulated 25 galactic-scale lenses (including 8 re-discoveries) so far and expect 
to eventually obtain a statistical sample of from 60-70 lenses@ These well-defined samples are 
particularly useful not only for constraining cosmological parameters such as the present-day 
matter density Om.n, dark energy dens ity Q Xj q and its equation of state w x (e.g. 



2002; Chae 2003 



Chae et al. 



Mitchell et al.ll2005l ) but also for const raining the s tatistical propert ies of 



galaxies such as optica l region velocity dispersions (e. 



galaxy evolutions (e.g. JChae fc Maoll2003 



Ofek et al. 



Chae 



20031 ). 



20051 : IChae et al.l 120061 ) and 



The sample from the completed CLASS, in particular its subsample of 13 lenses strictly 
satisfying w ell-defined selection criteria (the CL ASS statisti c al sam ple; iBrowne et al.ll2003l ; 
Chae 2003 ). was first extensively analyzed by Chae et al. ( 2002 ) and Chae ( 20031 ). who 
found fl m ,n ~0-3 ass uming a flat cosmology and adopting non-evolving galaxy populations. 



Mitchell et al.l (120051 ) re-analyzed the CLASS statistical sample based on the velocity disper- 
sion f unction (VDF) of early -type gala xies directly deriv ed from th e SPSS Data Release 1 
(D Rl:lstoughton et alJ liooa) galaxies Jsheth et al.ll2003h . However, Ichael (boosh finds that 
the Isheth et al.l (hooah VDF of early-type galaxies would imply a significantly underesti- 
mated abundance of early-type galaxies based on th e Wilkinson Microwave Anisotropy Probe 
(WMAP ) 1st year cosmolo gy ( Spergel et al .1120031 ) and the CLASS statistical sample. Just 
recently, IChoi et al.l (120071 ) have made a new measuremen t of the VDF of early-type galax- 



ies based on the much larger SDSS Data Release 5 (DR5; lAdelman-McCarthv et al. 



galax ies em ploying a new an d more reliable method of classifying galaxies (jPark fc Choi 



2007 



20051 ). The IChoi et al.l (120071 ) VDF has a much higher comoving number dens ity of early- 



type g alaxies and a different shape for the lower velocity part compared with th e lSheth et al. 
(120031 ) VDF. The lChoi etaD (l2007f ) early-type number density is in favor of the|Chae] (bood ) 
results. 

The goal of this work is to improve strong lensing statistics using the SDSS DR5 VDF 
of early-type galaxies. Our focus shall be to put independent constraints on f2 m and w x 
assuming a flat cosmology. We shall consider both no evolution and a evolution of galaxies 



1 lit tp :/ /cf a- www . harvard . edu/castles 

2 N. Inada, 2006 SDSS collaboration meeting in Seoul (http://astro.snu.ac.kr/ sdss/). 
3 The actually used data set is called DR4plus which is very similar to the DR5. 



- 3- 



based on the pre diction by a semi-analytical model of galaxy formation (IKang et al.l 12005 



Chae et al.ll2006l ). In §2, we briefly describe the data and the analysis method. We present 
and discuss the results in 83. 



Data and Method 



The co moving number density of galaxies as a function of velocity dispersion (a) can be 
descr ibed by the modified Schechter function 4>(a) given by (jSheth et al.ll2003l ; iMitchell et al. 
2005h 

(3 da 



dn 



{a)da = (f>* 



a 



cxp 



1) 



V{a/(3) a 

where 0* is the integrated number density of galaxies, <r* is the characteristic velocity dis- 
persio n, a is the low-velo city power-law index, and (3 is the high- velocity exponential cut-off 



index. Sheth et al. (the number density being updated by lMitchell et al.ll2005l ) found 

from the SDSS DR1 for the early-type galaxy population 



<?*, «, 0)dri = [(4.1 ±0.3) x 10~ 3 h 3 Mpc" 3 , 
88.8 ±17.7 kins" 1 , 
6.5 ±1.0, 1.93 ±0.22], 



(2) 



where h is the Hubble constant in units of 100 km s 1 Mpc 3 . This was the first direct 



measurement of the VDF of early-type galaxies. Just recently, IChoi et al.l (120071 ) have found 
from the much larger SDSS DR5 



a, (3) 



DR5 



[8.0 x 10~ 3 h 3 Mpc" 3 , 

161 ±5 km s~\ 

2.32 ±0.10, 2.67 ±0.07]. 



(3) 



(Sheth et al. 


2003; 


Choi et al. 


2007) 



in Figure Ufa). The DR5 VDF is clearly quite different from the DR1 VDF both in the 
number density and the shape for the lower part of a. This is in large part due to the 



improved galaxy classification scheme of iPark fc Choil (120051 ). who makes use of a SDSS u-r 
color versus g-i color gradient space. 

While early-type galaxies dominate strong lensing, late-type galaxies cannot be ne- 
glected. Among the radio -selected len ses with known galaxy types from the CLASS and the 
PANELS (see Table 1 of IChael 120051 ) . about 30% are late-types. For the late-type galaxy 
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population, the direct measurement of the VDF is complicated by the significant rotations of 
the disks. We proceed as follows. We ado pt the luminosity function (LF) of late-type galax- 
ies from the SDSS DR5 (IChoi et al.ll20071). We turn the LF into a circular velocity function 
using a Tully-Fisher relation (iTully &: Piercd l2000l ). Finally, we turn the circular velocity 
function into a VDF assuming that the circular velocity is proportional to the inner velocity 
dispersion as would be the case in an isothermal model. In principle, we can estimate all the 
parameters of equation ([TJ) for the late-type population using the Tully-Fisher relation and 
a galaxy model. However, we shall leave a* free and determine it from the image separation 
distribution. Our adopted VDF for the late-type galaxy population is then 



a, (3) 



late 



[1.13 x 1CT 1 h 3 Mpc" 



c4 late) , 0.3, 2.91]. 



(4) 



A late-type VDF with a 



(late) 



133 km s determined from lensing in §3 is shown in Fig- 



ure Ufa). It is interesting to n otice that this VDF of late- type galaxies matches relatively 



well that of ISheth et all (j2003l ) who determined using a Tully-Fisher relation taking 



into account the scatter of the Tully-Fisher relation. 

The SDSS galaxy populations refer to a redshift range of < z < 0.2 while the radio- 
selected lenses are in the range of 0.3 < z < 1, so that galaxy evolutions must be taken into 
account. Most of previous works have been done with the assumption of no ev olution of ea rly- 
type galaxies from z — 0-1 relying on several observational arguments (see IChadl2003l and 
reference s therein) such as the fundamental plane and the color-magnitude relations. More 
recently, iMitchell et all (120051 ) have taken into consideration galaxy evolutions using the 
prediction by the exten ded Press-Schechter mo del of structure formation which is calibrated 
by N-body simulations (ISheth fc Tormenll2002l ). We take into consideration galaxy evolutions 



using the predict ions by a recent semi-analytical model of galaxy formation (IKang et al. 



2005 



Chae et al.ll2006l ) which is based on the high- resolution N-body simulations of iJing &: Suto 
( 120021 ). Specifically, we constrain t he evolutions of th e VDFs using the evolutions of the virial 
circular velocity functions used by lChae et all (120061 ) assuming a power-law relati on between 



the vi rial circular velocity and the velocity dispersion as given by equation (6) of IChae et al. 
( 120061 ). We consider the evolutions of both the number density 0* and the characteristic 
velocity dispersion a* as follows: 



J,(z) = 0* )O (1 + z) Vn \ a*{z) = o-*,o(l + zf\ 



(5) 



We obtain th e following best-fi t parameters for the early-type and the late-type galaxy 
populations of IKang et all 12005 : 



(-0.229, -0.01) for early-type and (1.24, -0.186) for late-type. 



(6) 
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We us e all the lensing properties of the CLASS and the PANELS lenses (Table 1 of 
Chad 120051 ) . Specifically, we use both the image separation distribution of all 'single' lenses 
as well as the number statistics of the well-defined CLASS statistical sample of 13 lenses. 
W e use the sam e singu l ar iso therm al ellipsoid lens m odel and analysis method recently used 
by I Chad (120031 ). IChad (120051 ). and IChae et al. (2006). The li kelihood function is essentially 
the same as that given by equation (11) of IChae et al.l (120061 ). namely 



(7) 



where Spis{j) is the relative image separation probability, p(k) is the total multiple-imaging 
probability due to both the early-type and the late-type populations, and 5p{l) is the dif- 
ferential lensing probability of the specific separation, lensing galaxy type, lens and source 
redshifts. The only difference here is that for the lenses of unknown galaxy morphologies the 
image separation probability Spis(j) is calculated as a weighted sum due to the early-type 
and the late-type galaxy populations where the weighting factors are the diffe rential lensing 
probabilities for the observed specific image separations (eq. 29 of IChad 120031 ). 



3. Results and Discussion 

We first consider a flat cosmology with an Einstein cosmological constant. Figure [0(b) 
shows the behavior of the likelihood function (eq. [7]) as the matter density Q m ,o is varied. 
The results shown in Figure [T](b) are based on the SDSS DR1/DR5 early-type VDF (eqs. [2] 
&ED and the late-type VDF (eq. H}. For the DR5 VDF we find ft m ,o = 0.25±g;J| and 0.26±° f 8 



70 CL) respectively for the evolving (eq. [6]) and non-evolving populations of galaxies. The 
fitted value of for the late-type population is 133 km s" 1 and 127 km s _1 respectively 
for the evolving and non-evolving cases. 

For the DR1 VDF we find ft m ,o = 0.18t^ 6 and 0.19i^ 6 (68% CL) respectively for the 
evolving (eq. [6]) and non-evolving populations of galaxies. The difference between the DR1 
and the DR5 results can be understood from the behavior of 0<r 4 , which is proportional to 
the differential lensing probability , shown in Figu r e [Ha) . Notice that our results based on the 



DR1 VDF are different from the iMitchell et al.l (120051 ) results because of differences in the 



calculations of magnification biases and cross sections satisfying observational selection func- 
tions. For example, for the differential number- flux density relation of \dN/d S\ cx (S/Sq)~ v 



with S = 30 mJy, we use r] = 1.97 (2.07) for S < S (S > S ) while Mitchell et al. 
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( 120051 ) uses erroneously r] — 2.1 for any Taking into account the flux ratio limit for the 
doubly-imaged systems in the CLASS statis tical sample, namely that the fa in ter-to-brighter 



i mage flux ratio must be greater than 0.1 (IBrowne et al.ll2003l ; IChaei 120031 ). iMitchell et al. 
( 2005h finds B = 3. 97 using the singular isothermal sphere (SIS) model, where B (eq. 15 of 



Mitchell et al.ll2005l ) is the magnification bias times the cross section satisfying the flux ratio 
limit divided by the unbiased cross section. However, we find B pa 3.36 for the SIS taking 
corr ectly into account th e CLASS observational selection functions. Another difference is 
that IMitchell et al.l (120051 ) uses the SIS while we use the singular isothermal ellipsoid (SIE). 
For example, a lens ellipticity of 0.4 can amount to a difference of Af2 mj o ~ —0.05 compared 
with the spherical case because of the variation of the magni fication bias and cro ss section 
for the equal numbers of oblates and prolates (see, however, iHuterer et al.l 120051 for other 
possibilities). 

The directly measured SDSS DR5 VDF of early-t ype galaxies is more relia ble than VDFs 
inferred fro m early-type LFs using a Faber- Jackson (IFaber fc Jackson! 119761) relation (e.g., 



Chadl2003l ) because of significant scatters in the relation ( Sheth et al 



VDF is also likely to be more reliable than the DR1 VDF (iSheth et al 



2003). The SDSS DR5 



20031 1 not only because 



the DR5 VDF is based on a much larger volu me sample but also b ecause it is based on a 
more reliable galaxy classification technique by iPark fc Choi! (120051 ). Therefore, the SDSS 
DR5 VDF in conjunc tion with galaxy e volution models f rom a recent semi-analytic model 
of galaxy formation (IKang et al.l 120051 : IChae et al.l 120061 ) removes in large part potential 
systematic errors in one of the main ingredients of strong lensing statistics. 

Our derived values of the matter density based on the SDSS DR5 VDF are in excellent 
agreement with the r esults from the WMAP and the large-scale structures in the SPS S 
luminous red galaxies (jSpergel et al.ll2003l . 120061 ; iTegmark et al.ll2004l ; lEisenstein et al.ll2005l ) . 
The SDSS lens search is eventually expected to discover from 60-70 strongly lensed quasars. 
Thus, we expect that the precision of strong lensing statistics will improve by a factor of two 
or better in the near future (in particular if the CLASS statistical sample is combined with the 
SDSS sample). Furthermore, with in a few decades next generation observa tion tools such as 



the Square Kilometre Array (e.g.. lBlake et al.ll2004l ; iKoopmans et al.ll2004l ) will improve the 



precision of lensing statistics by orders of magnitude. Perhaps, strong lensing statistics may 
play even more important roles in the future for uncovering the physical processes of galaxy 



formations and evolut i ons (see, e.g.. iKochanek fc White! 12001 



2003; Ofek et al. 2003: Chae et al. 2006). 



Keeton 2001: Chae fc Mao 



4 While the lensed sources have observed flux densities S x> > 30 mJy, the magnification factor is given 
as an integration for the ra nge from S = S b to S h/[J-m&x where /i ma x is the maximum possible theoretical 
magnification (see eq. 36 of lChael l 2003). Thus, it is essential to use rj = 1.97 for S < 30 mJy. 
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Finally we consider a flat cosmology with a generalized dark energy Q x> q with a non- 
evolving equation of state w x . Figure [2] shows the confidence limits in the fl m ,o-w x plane. We 
have w x < —1.2 at the 68% CL (w x < —0.5 at the 95% CL), so that strong lensin g appears 
to marginally favor a super-neg ative equation of st ate (e.g., ICaldwell et al.ll2003l ). On the 
other hand, the WMAP results (ISpergel et al.ll2006l ) favor w x = — 1 (Einstein's cosmological 
constant). It is not clear at the present why lensing data appear to favor w x < — 1. It will 
be interesting to see whether this remains so with future lensing data. 
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Fig. 1. — (a) Thick lines are the DR5 (solid) and the DR1 (dashed) early- type VDFs and the 
late- type VDF (dotted). Thin lines show the behavior of 0<r 4 (which is proportional to the 
differential lensing probability) for these VDFs. (b) The behavior of the likelihood function 
(eq. [7]) as fl m ,o is varied for a flat cosmology with an Einstein cosmological constant. The 
solid (dashed) lines are based on the SDSS DR5 (DR1) VDF of early-type galaxies and the 
SDSS DR5 LF of late-type galaxies. The black and red curves correspond respectively to 
the cases of assuming no evolution of galaxies and the evolutions of galaxies predicted by a 
recent semi-analytical model of galaxy formation. 



- 11 - 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

CI n (= 1 - G n ) 
m,U v x, U' 

Fig. 2. — Constraints in the plane spanned by the matter density (= 1 — &>x,o in our 
assumed flat cosmology) and the non-evolving equation of state w x of a generalized dark 
energy based on the SDSS DR5 VDF of early-type galaxies and LF of late-type galaxies 
through radio-selected strong lensing statistics. We have assumed the evolutions of galaxies 
predicted by a semi-analytical model of galaxy formation. 



